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INTRODUCTION 


In  the  first  volume  of  these  Reports*  a  discussion  and  tables  have  been  given  describing  the  properties 
of  a  supersonic  flow  around  cones  when  the  direction  of  the  stream  is  exactly  parallel  to  the  axis  of  the  solid. 
Such  conditions  can  be  realised  experimentally,  for  instance,  in  a  wind  tunnel  where  a  model  exposed  to  the 
stream  car,  be  rigidly  maintained  in  •  fixes]  position,  but  not  in  the  much  more  important  practical  case  of  a 
solid  travelling  at  supersonic  speed  through  air  at  rest:  for  any  actual  projectile  in  free  flight  wilt  necessarily 
yaw,  and  the  yawing  motion  superposed  upon  ordinary  translation  will  bring  about  deviations  from  n  purely 
axial  flow,  the  amount  of  which  must  be  investigated.  If  the  missile  is  reasonably  stable,  its  yaw  may  be 
expected  to  remain  small;  and  the  possibility  of  modifying  an  axial  flow  by  appropriate  perturbation  method* 
to  account  for  the  effects  of  a  small  yaw  suggests  itself. 

An  approximate  method  along  these  fines,  Past  J  on  thr  assumption  of  isentropy  and  irrotationaiity  of 
the  stream  around  a  slightly  yawing  cone  has  been  developed  by  Tsien,*  and  a  more  refined  method,  based, 
however,  on  similar  simplifying  assumptions,  has  later  been  given  by  Sauer.*  With  the  outbreak  of  the  recent 
war  the  interest  in  the  flow  of  air  around  yawing  conical-headed  projectiles  became  at  once  acute  and  stimu¬ 
lated  renewed  attacks  at  more  general  types  of  solution  of  the  ur..'»rlying  hydrodynamical  problem.  As  « 
result,  an  investigation  based  on  less  restrictive  simplifying  assumptions  has  been  produced  by  Kt:ush  and 
Critchfield,*  and  their  method  was  later  greatly  perfected  by  Stone.*  In  January'  1943,  ar  advance  version  of 
Stone’s  work  was  brought  to  the  writer’s  attention  by  Dr.  Raymond  J.  Seeger,  then  of  Section  Reie  cf  the 
Bureau  of  Orduanee  of  l.’.S.  Navy  Department,  with  the  request  that  its  consequences  be  investigated  quanti¬ 
tatively.  The  project  was  subsequently  approved  by  Commander  Charles  C.  Bramble,  then  acting  as  the 
liaison  officer  between  the  U.S.  Naval  Proving  Ground  at  Dahlgren,  Virginia,  and  our  group.  A  substantial 
part  of  the  resuits  obtained  in  the  course  of  this  investigation  in  the  years  1943-  4?  is  now  being  presented  in  this 
volume.  Since  the  hydrodynamical  theory  which  was  used  as  a  basis  for  our  computations  has,  however,  not 
yet  been  made  public,  it  seems  appropriate  to  open  this  book  with  a  brief  outline  of  the  theory  and  a  discus¬ 
sion  of  its  merits  and  limitations.  In  what  follows  such  an  outline  will  be  given. 

EQUATIONS  OF  THE  PROBLEM 

Consider  (cf.  Fig.  1)  a  rectanp-tlar  system  of  coordinates  (1,  2,  3),  where  r,  0  denote  the  usual  spherical 
polar*.  The  equation  of  a  circular  cone  with  the  vertex  at  the  origin,  and  the  axis  making  an  angle  <  with  the 
coordinate  axis  (1)  can  be  expressed  as 

r  cos  6  cos  *  r  sin  9  sin  «  cos  0  =  r  cos  9,  (1) 

where  9,  denotes  the  semi-apex  angle  of  the  solid  cone.  Let  us,  in  what  follows,  consider  the  angle  of  yaw  < 
to  be  a  quantity  so  small  that  its  squares  and  all  higher  powers  can  be  ignored  To  this  order  of  accuracy  the 
equation  (1)  of  the  cone  can  be  reduced  to 

9  —  6,  +  *  cos  0  +  -  -  -  -  ,  (1.1) 

and  our  problem  will  be  to  specify  the  profierlies  of  a  supersonic  air  stream  around  such  a  cone  to  the  same 
degree  of  approximation. 

Let,  as  usual,  p  and  p  denote  the  pressure  and  density  of  the  gas  at  any  point  of  the  flow.  The  equations 

\.  Super  eunv:  Flow  of  Air  around  Come,  Massachusetts  Institute  (4  Technology,  Cambridge,  Mtw  ,  1M7. 

i  Joum.  Aero.  Sci..  5,  480,  1058. 

3  Zeitjchr.  f.  Luf  tfohrtf orach  ung.  1*.  148,  1644. 

*  The  Preteure  on  a  Cone  moeine  with  Small  Fair  ai  High  I’dotUu,  NDRC  Report  Div.  A,  1644. 

5.  The  Aerodynamics  of  a  Slightly  Yawimf  Supereonic  (  one,  NDRC  Report  Div.  1,  164*.  'Hie  writer  understand*  from  a  recent 
letter  by  Dr  Stone  that  t  full  dreflo  version  of  hi*  theory  may  no on  be  published  in  •  wientifir  pe  nodical 
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of  motion  ol  a  perfect  fluid  in  a  steady  state  then  assume  the  form 

1*E  _ 

p  dr 
]  dp 
pr  d6 
1 


du  ,  *  9u  t  u>  du 

u —  -f - -4 - — - h 

dr  r  d6  r  »n  6  d$ 


dr 
u  — 
'  dr 


dir  r  dir 

u - (-  —  ■ 

dr  r  56 


dv 


r»+tJ 


=  0. 


r  dr  it  dr  ^ 

r  56  r  sin  6  dip 


»r  —  v1  cut  6 

+  - ...  . ---  —  0  , 

r 


H - : — 


r  sin  8  54>  pr  sin  8  dip 


dp  .  tns-t-rv  cot  8 
—  +  -  ■  =0, 


I 

I 


(*) 


white  the  equation  of  continuity  can  be  expressed  as 

Add 

—  (prhi  sin  it  +  —  (per  sin  6)  -t - (pnr)  —  0 ,  (8) 

dr  3ff  tty 


where  u,  e,  w  denote  the  velocity  components  in  the  direction  of  increasing  r,  8,  and  0.  Let  us  seek  now  such 
solutions  of  these  equations  which  are  independent  of  r,  end  expand  »,  r,  ir  as  well  as  p  and  p  in  series  of  the  form* 


n  =  si  -f-  «  £  xm  cos  rty  +  -  .  .  . 

«-0 


r  «  f  cos  nj>  4  -  •  •  • 

«-o 

tr  =  c  2  tm  sin  ti4>  4  ‘  -  *  • 

fH 

« 

p  =  i  +  <  J  t,  cos  n0  +  •  *  ■  • 

II —0 


P  =  p  4  «  2  t*  cos  4  •  ’  *  * 

A  -0 


(4) 


where  6,  f,  p,  and  p  stand  for  the  respective  quantities  characteristic  of  the  purely  axial  flow,  and  x,  y,  z,  g, 
and  £  specify  the  first-order  effects  of  yaw.  The  velocity  components  H  and  f  are  known  U,  satisfy  the 
following  differentia]  equation 

<Pu  .  .  a’(«  +  t  cot  9) 

—  +  M  = - - . -  (S) 


of  Taylor  and  Maccoll,’  where 


and 


(«) 


—  l)(c=—  *i?  —  r») 

2 


denotes  the  square  of  the  local  velocity  of  sound;  y.  the  ratio  of  specific  heats  of  the  streaming  gas;  and  r,  the 
velocity  which  the  gas  would  attain  if  allowed  to  expand  adiabatically  into  a  vacuum.  As  ali  barred  quantities, 
also  z,  y,  i,  ij,  and  £  will  similarly  hereafter  be  supposed  to  be  functions  of  8  only,  and  their  forms  sought  as  solu¬ 
tions  of  ordinary  differential  equations  with  8  as  the  independent  variable.  Inserting  (4)  in  (2)  and  (8)  we  find 

8.  Justification  lor  neglecting  the  line  terms  in  tbene  Fourier  expansions  will  be  found  in  Stone's  forthcoming. paper  oo  the  subject 
It  should  softer  here  to  #ay  that  such  lermt  could  not  meet  the  appropriate  boundary  rendition*  on  the  rone. 

7.  Proc.  Roy.  Soc.,  A,  1 M,  478,  1®S3.  Cf.  also  the  introduction  to  volume  1  of  these  Report*. 


IX 


indeed  tiie  m  w  f  auctions  to ! 


deSfifu  Ly 

V.  -=*•', 

ty.'4-(*"+*)».+*r.+*,Vp-^V^-o. 

9*»'4-  (il+f  cot  •)*,— ne./p  sin  5  =0, 

(So*  g)'}  nr.  4-  4r.  +  nr.  cw  »  +  S(t.'g)' «0, 


!  cot 


(7) 

(8) 
(8) 

(10) 


w  litre  hcCcuu  denote  diffeieatiatioiu  with  respect  tc  !. 

The  equatior  of  stair  relating  p  and  p  remains  to  be  pfeified.  If  there  were  do  yaw,  the  ratio  p/py 
should  be  constant  everywhere  behind  the  shock  wave  (if  any),  though  different  from  that  in  front  of  it.  The 
appearance  of  yaw  will  change  inis  situation,  rince  a  shock  wave  would  alter  the  entropy  by  amounts  which 
vary  with  the  angle  of  incidence  and  therefore  with  <t>  If.  however,  we  agree  to  ignore  the  effects  of  viscosity 
and  best  conduction  in  the  Hot  behind  the  shock  wave,  it  can  be  proved*  that  the  flow  will  remain  adiabatic  at 
local  for  each  individual  at.  particle  Since  the  motion  is  assumed  to  be  steady,  the  actual  paths  of  the  particles 
coincide  with  the  mathematical  streamlines:  therefore,  for  small  yaw,  the  equation  of  state  along  each  stream¬ 
line  should  semribh  reduce  to 


(11) 


1  dp  y  dp 

p  at  ~  P  ae * 

which  is  integrate  with  respect  to  8.  It,  ultimately,  we  insert  in  (11)  the  above  expansions  for  p  and  o,  we 
find  that,  to  the  first  order  ip  small  quantities  the  equation  of  state  can  be  reduced  to 


(18) 


where  dm  is  a  constant.  Equations  (7) —(10)  and  (14)  together  with  their  appropriate  boundary  conditions 
contain  a  complete  specification  of  our  problem. 

The  boundary  conditions  are,  in  general,  the  same  as  for  a  purely  axial  fiow :  namely,  at  the  surface  of  the 
solid  cone  the  normal  velocity  component  must  vanish,  whereas  at  the  shock  wave  the  standard  Rankine- 
Hugooiot  conditions  must  be  met.  It  can  be  shown  that  a. requirement  of  disappearance  of  the  normal 
velocity  component  at  is  tantamount  to  a  statement  that,  if 


1,  pi  -  48, 


while  if 


«  m  l,  y.=0. 


(IS) 


It  follows  that,  if  equations  (4)  are  to  satisfy  this  requirement,  all  variables  with  subscripts  other  than  1  arc 
identically  iero.  In  consequence,  x'l  subscripts  will  hereafter  be  discarded.  To  the  degree  of  accuracy  we  arc 
worldug,  the  equation  of  the  shock  wave  then  assumes  the  form 


#  =  >»  +  »  cos  d. 


(14) 


where  i,  is  the  undisturbed  shock-wave  angle,  and  A  is  an  undetermined  constant.  Inserting  this  in  the  Ranlrine- 
Hugoniot  shock-wave  conditions  (guaranteeing  the  continuity  of  the  tangential  velocity  component,  of  energy, 
and  of  impulse)  we  find  that,  at 

(15) 

(1«) 

(17) 

(18) 
(19) 


A  =  —  z/(B  -|-  U  sin  9), 
x  +  *  sin  t  =  0, 
ix  cot  8  -r  y  +  8i/p  =  0, 

—x  cot  f +  >  +  1| /pi  =  0, 

*  {»(log  p)‘  +  *9/<_y  ■+•  1)}  +  y(t  +  i1  tan  0)  =  0, 

where  V  denotes  the  velocity  of  the  undisturbed  stream  n  front  of  the  shock  wave. 
8.  Cf„  frtr  lmUidcc,  V Maonji,  Quart.  of  Appl.  M-»th  J,  1$,  1*45. 
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SOLUTION  OF  THE  EQUATIONS 

Equations  (15)  (19)  together  with  (15)  specify  the  requisite  boundary  conditions  of  our  problem.  In 
order  to  proceed  with  its  solution  let  u#  eliminate,  in  equation  (8),  y  and  {  by  mefcos  of  (7)  and  (14) ;  we  obtain 

ax"4-(«"4-  «)r'f  fa  +  --  £-h-d\~G.  (40) 

P  7P  <r  > 

This  equation  can  be  integrated  once  and  yields 

ftr'  Hr  4-  —  -| — — -4-fc~0,  (41) 

p  r-i  i 

where  k  is  a  constant.  If  we  set,  in  this  equation,  0=*0„  an  insertion  of  (17),  (18),  and  (14)  into  (41)  proves 
that  1  =  0.  Hence, 

#*'+  Uz-h-+~  &  “0,  (44) 

p  7 — 1  P 

which  represents  an  analogon  to  the  Bernoulli  equation,  and  could  in  fact  have  been  derived  from  it.* 

Neat  combine  (9)  and  (44)  to  obtain  a  differential  equation  for  x+e  sin  8  of  the  form 

®(x  +  a  sin  8) '  4-  U(x  4-  *  sin  ®)  4 - —  ^  *  0 .  (4S) 

y-1  i 

This  is  a  orst-order  linear  equation  whose  integrating  factor  is  (— Pp  sir.  9)”*;  its  solution  takes  the  form 

.  <fV~fp  sin  8  f*  _ pde 

y  —  1  J ’«  —  tp'i ~  ip  sttl  <? 

where  the  limits  of  integration  have  been  adopted  so  as  to  make,  in  accordance  with  (18),  x-f-i  sin  8  vanish  on 
the  shock  wave.  Now  equations  (?),  (12),  (22).  and  (24)  enable  all  unknowns  other  than  x  to  be  eliminated 
from  (10).  Going  through  the  requisite  and  rather  troublesome  algebra  we  can  ultimately  rewrite  (10)  in  the 


x"/d  +  A(x'/d)  +  B(x/d)  +  6*0, 

A  *»  cot  8  +  x[2u  +  S®  cot  8  -4*  (7  4*  l)Xf(i<  4"  ®  cot  9)1 , 
B  «  1  -  cot*  8  +  x[x(y  -  l)fi(tf  4-  0  cot  8)  -  ®  cot1  e] , 

„  1  +  W  , — -T-r  (•  1* 

(t  1)  sin  «  Jtw  -fpV  _  it:  sin  $ 


subject  to  the  boundary  conditions  that,  when  0=fl» 


y(y~l)(“  4-  scot®) 


x'  _  2a*  cot  8  (2a  —  ®  cot  9) 

d  7(7*—  !)(■**  4-  foots)’’ 

whereas  if  8—8,, 

x‘  -  m,. 

B.  Out  being  tsken,  however,  to  «ilow  for  the  feci.  th»t  the  present  flow  is  neither  irroutiooaj  nor  isentropio. 


TABLES  OF  SUPERSONIC  FLOW  AROUND  YAWING  COSES 

Once  equation  (25)  with  the  above  boundary  condition*  ha*  been  integrated  —  numerically  or  otherwise 
—  the  solution  of  our  problem  ia  complete,  for  all  other  unknown*  required  for  a  uniqtx  description  of  the  flow 
can  be  obtained  from  i  or  it*  derivative*  by  a  aerie*  of  purely  algebraic  operation*.  It  transpire*  that  the  thock 
wave  attached  to  a  slightly  yawing  cone  in  supersonic  flight  will  continue  to  be  a  circular  cone,  of  the  same  apex 
angle  a*  in  the  non-vaw  case,  but  with  a  yaw  i  of  it*  own,  given  by 

i  —  at, 

where 

a  — _ a'd  Cot  * -  ~  £  .  (»,,» 

7(7-  !><*»  +  Scot#)1  4R-»eot* 


evaluated  at  0  »  0..  Furthermore,  it  is  evident  that  the  plane  of  yaw  of  the  shock  wave  will  be  the  same  as 
the  plane  of  yaw  of  the  cone, 

The  forces  acting  on  a  slightly  yawing  cone  are  ail  expressible  in  terms  of  *.  In  particular,  the  pressure 
at  a  point  P(r,  6,  +  •  cos  <p,  pi  of  the  cone  can  be  written  as 

p-A  +  t.  «cosp.  (28) 

where  p  (ef.  Fig.  1)  is  the  aaimuthai  angle  measured  from  the  plane  of  yaw  around  the  axis  of  the  cone.  The 
norma!  force  Ny  on  •  cone  of  finite  rise,  and  the  distance  ktl  from  the  vertex  to  the  center  of  head  pressure 
follow  readily  by  integration:  we  obtain 


and 


(SO) 

ha  «•  —  l  sec1 1, , 

3 

(31) 

respectively,  where  D  denotes  the  diameter  of  the  base  of  the  cone  and  i,  the  distance  from  the  vertex  to  the 
base.  The  coefficient  K#  of  the  force  acting  on  the  cone  in  the  direction  perpendicular  to  its  motion  ia  then 
defined  by 


-V* 


„  I  (54) 

pJPIPm  *  8  mU> 

while  the  coefficient  KD  of  the  head  drag  acting  on  the  cone  in  the  direction  of  its  axis  is  found  to  have  a  value 


Kd  “  J  (A  -  ?.)£*.  (33) 

which  is  identical  with  that  fora  purely  axial  flow.  In  these  equations  pt  and  />>  denote  the  pressure  and 
density  of  undisturbed  air  in  front  of  the  shock  wave. 

It  may  be  added  parenthetically  that  the  above  equations  (28)  and  (34)  should  continue  to  hold  good  as 
long  at  terms  of  the  order  of  «>  and  higher  remain  negligible.  When  squares  of  the  angle  of  yaw  are  taken  into 
account,1*  there  arises  a  contribution,  du*  to  yaw,  to  the  head  drag,  and  the  shuck  wave  cease*  to  be  a  circular 
cone  (ita  ha  sis  being  defined  by  a  curve  of  the  sixth  degree) :  but  the  yaw  of  the  shock  wave  and  the  normal  force 
coefficient  remain  unaffected. 


DESCRIPTION  OF  THE  TABLES 

Equation  (25)  given  above  represents  the  fundamental  equation  of  our  problem  in  much  the  same  way  as 
equation  (5)  has  been  our  basis  for  investigations  of  axial  supersonic  flow.  It  is  a  second-order  nonhomogeneous 
differential  equation  with  variable  coefficients,  the  solution  of  which  cannot  obviously  be  expected  to  be 
expressible  in  a  finite  number  of  terms.  The  problem  for  computation  therefore  consists  in  solving  (25) 

10.  These  ejects  have  also  beer  invest  igsted  by  Dr.  9tme  sad  srr  oontBiard  in  •  paper,  now  in  memucript,  which  the  present  writer 
wee  pMqri  to  see  ia  advance  of  piABeethm. 
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numerically  (or  x/d,  starting  it  J-l,  with  the  initial  condition*  (88)  and  intagrating  backward*  until  the 
•olid  cone  U  reached ;  the  inner  boundary  condition  (87)  then  determine*  the  value  of  4.  Once  thia  hr*  bean 
accomplished  and  x  with  x’  tabulated  a*  (unction*  of  4,  the  remaining  unknowns  ),  *,  q,  and  t  can,  aa  already 
mentioned,  be  expressed  in  terms  of  x  and  x'  by  a  series  of  purely  algebraic  operations.  In  more  explicit  terms, 
equation  (7)  asserts  that  y  — x'  and,  by  (84), 


£  V  tin  f  _  -r 

d  1  4-  Xi  dm in  #  ’ 

while  (remembering  that  jf/p  =  a’/r)  the  ‘'Bernoulli  integral”  (88)  yields 

1  m  —  ~  —  T(**  +  **') 

f  y-l  o> 

and,  from  the  equation  of  state  ()8)  it  follows  that 


(«) 

(M) 


p 


In  order  to  be  able  to  integrate  equation  (85)  numerically,  we  first  have  to  evaluate  its  coefficients  A,  B 
and  C  a*  functions  of  9.  These  coefficients  are,  in  turn,  known  functions  of  the  velocity  components  4  and  9 
of  the  corresponding  non-ysw  case,  the  values  of  which  were  published  extensively  in  Volume  1  of  these  Reports, 
The  coefficients  A  and  B  are  purely  algebraic  functions  of  t?  and  f ,  while  C  involves  an  additional  quadrature." 
Tables  of  all  three  coefficients  have  been  prepared  on  the  basis  of  the  velocity  components  4  and  9  which  them¬ 
selves  were  computed  to  six  decimal  figures,  and  the  integration*  of  (85)  for  x/d  carried  out  to  five  figures.  The 
remaining  unknowns  y,  i,  *,  and  £  were  then  computed  from  equations  (7),  (34),  (35),  and  (18),  and  eventually 
rounded  off  (in  most  cases)  to  four  decimal  figures  which  are  regarded  to  be  reliable  within  one  unit  of  the  last 
place.  More  decimals  were  retained  only  in  certain  solutions  for  slender  cones  where,  in  the  majority  of 
values  given  in  some  columns,  the  decimal  point  waa  followed  by  one  or  more  sera*. 

Fart  I  of  this  volume,  which  follows  the  introduction,  contains  the  results  of  180  such  integration*, 
carried  out  for  a  value  of  7^  1.405  which  appears  to  be  the  best  existing  approximation  to  the  ratio  of  specific 
heats  of  dry  air.  It  includes  18  groups  of  integrations  for  cones  of  semi-apex  angles  ranging  from  5*  to  50*; 
and  within  each  group  4,/c  increases  in  regular  increments  from  0.3,  which  is  well  within  the  domain  of  the 
“second"’  solutions  of  equation  (5)"  (marked  with  **S"  following  the  particular  value  of  the  Mach  number 
U/ai  in  the  heading  of  a  tabic),  up  to  the  limit  for  which  the  Mach  number  becomes  infinite.  Successive 
columns  of  the  individual  tables  contain,  respectively,  the  quantities  x,  y ,  *,  q/f,  and  f/g,  by  means  of  which  the 
actual  velocity  components  u,  v,  w,  the  pressure  p.  and  the  density  p  at  any  point  between  a  slightly  yawing  cone 
in  supersonic  flight  and  its  shock  wave  can  be  found  from  equations  (4) .  As  the  velocity  components  u  and  9  in 
Volume  1  of  these  Reports,  also  the  present  quantitu *  x,  y,  and  t  halt  been  expressed  in  lermt  of  the  velocity  c 
which  the  air  would  attain  if  allowed  to  expand  adiaba/icaUy  into  a  vacuum.  The  absolute  value  of  c,  for  a  given 
stream,  can  easily  be  found  from  ti>e  equation 

«*-P»  +  —  .  (88) 

7- 1 

where  U  denotes,  as  before,  the  velocity  of  the  undisturbed  stream  (equal  to  the  actual  speed  of  a  cone  through 
air)  and  alt  the  velocity  of  sound  in  the  undisturbed  stream  which,  for  perfect  gas,  is  likewise  a  known  function 
of  absolute  temperature.  T'he  absolute  '-slues  of  x,  y,  or  1  at  any  point  behind  the  shock  wave  may  then  be 
found  by  multiplying  the  tabulated  non-dimensional  quantities  by  the  ahsolute  value  of  r  following  from  the 
above  equation  (38)  and  appropriate  to  the  given  particular  case. 

1 1 .  A  check  00  this  quadrature  is  aCrwtiad  by  the  requirement  (based  on  equation  |*S|)  that,  on  tbs  cone, 

c-  _  —*(*)  - 

L  y-l\pij,  tycilin'V 

It.  For  ft  rHrwino  of  the  meffioinf  of  tike  tod  “aftooevi”  aototioo*  <f.  i be  introduction  to  Voter**  1  of  tteae  Report*. 
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A  selection  o f  the  interval  of  tabulation  for  all  solutions  presented  in  Part  I  has  been  guided  by  the  deeire 
to  suable  the  user  to  interpolate  9-wi*e  for  intermediate  entriea  of  any  tabulated  quantity  with  the  aid  of  pro* 
portioeal  part*  alooe  In  order  to  make  thia  poaaible,  the  *eoond  difference*  of  all  tabulated  quantities  have, 
m funeral,  been  kept  under  20  unit*  of  the  U*t  place  by  a  suitable  choice  of  interval.  In  conaequenee,  the  error* 
committed  by  a  linear  interpolation  between  any  two  neighboring  entries  should  not.  in  general,  exceed  one  unit 
of  the  tart  tabulated  place.’*  The  interval  ct  tabulation  as  shown  in  the  tables  does  not,  therefore,  necessarily 
correspond  to  the  one  in  which  the  integration  was  originally  carried  out.  In  preparing  the  tables  for  publication 
the  columns  of  ail  quantities  were,  after  their  evaluation,  differenced  mechanically  using  two  National  Cash 
machine*  by  a  technique  described  by  Connie”  as  a  safeguard  against  any  accidental  error.  After  all  obvious 
slip*  had  thus  been  eliminated,  the  entries  were  punched  on  card*  (the  tine  for  each  9  on  a  separate  card),  tabu¬ 
lated.  and  proofread  with  greatest  care  against  the  original  work  sheet*.  The  corrected  set s  of  cards  were  then 
printed  by  an  I  B  M.  tabulator  on  prepared  forms,  and  tables  eventually  reproduced  by  photo-offset. 

The  data  presented  in  Fart  I  permit  one  to  determine  the  velocity,  pressure  and  density  of  air  at  any 
point  between  a  slightly  yawing  cone  in  supersonic  flight  and  its  shock  wave.  A  comprehensive  summary  of  the 
principal  physical  characteristics  of  these  solutions  is  given  in  Part  II,  the  successive  columns  of  which  indicate; 
(1),  the  adopted  velocity  component  if,  along  the  solid  surface  and  expressed  in  terms  of  <•;  (2)  and  (S),  the 
corresponding  semi-apex  angle  9,  of  the  shock  wave  and  the  Mach  number  17/«i“  Af ;  (4).  the  value  of  d  occur¬ 
ring  in  the  equation  of  state  (It)  and  specifying,  incidentally,  the  amount  of  entropy  change  across  the  stream* 
Hats;  (A),  the  shock-to-oone  yaw  ratio  l/«  as  defined  by  equation  (28.1);  (8),  the  coefficient  of  bead  drag  KD 
aa  defined  by  equation  (58);  and' (7).  the  coefficient  of  normal  drag  Km  following  from  equation  (92).  The 
values  given  in  columns  (1}-(S)  are  identical  with  those  already  published  in  Part  II  of  Volume  1  of  these 
Report*  —  save  for  those  few  tines  where  the  values  of  u,  in  Volume  1  are  prefixed  with  an  asterisk.  The 
respective  integration  of  axial  flow,  carried  out  to  only  five  decimal  figures  for  Volume  i,  was  subsequently 
drat*  over  to  six  decimal*  for  the  purpose  of  yaw,  and  the  outcome  indicated  sometime*  a  slight  change  in 
the  formerly  published  values  of  9U  and  (or)  M.  As  the  reader  can  verify  by  direct  comparison,  such  changes 
termed  out  to  be  barely  significant  (one  or  two  units  of  the  last  place)  —  which  reaffirms  our  confidence  in  die 
statements,  made  in  the  introduction  to  Volume  1  of  these  Reports,  concerning  the  general  accuracy  of  our  work. 
Column  (8),  i*»ntoiw«»g  the  values  of  head  drag  A'*>,  calls  for  an  additional  word  of  explanation.  This  quantity 
has  also  been  published  before  in  Part  XI  of  Volume  !.  For  the  present  work  it  was  multiplied  by  s  normalising 
factor  of  v/4,  in  accordance  with  equation  (88),  so  as  to  make  it  consistent  with  the  genera]  ballistic  usage  and 
with  values  of  Km  which  have  been  normalised  in  the  same  way. 

As  ha*  already  been  mentioned,  a  small  yaw  of  a  solid  cone  moving  at  supersonic  speed  will,  to  the  first 
t>-der  h\ small  quantities,  make  the  cone  of  the  shock  wave  yaw  at  an  angle  i  which  is  generally  different  from  «, 
and  give  nae  to  a  force  acting  on  the  solid  in  the  direction  perpendicular  to  its  motion  The  subsequent  Parts  III 
and  IV  of  the  present  volume  contain  detailed  tabulations  of  6/  <  and  of  Km  in  terms  of  9,  and  H$/c,  while  the 
diagrams  in  the  concluding  Part  V  represent  their  variation  graphically.  The  reader  who  wishes  to  look  them 
up  for  a  particular  pair  of  9,  and  the  Mach  number  M  rather  than  may  do  so  with  the  aid  of  a  conversion 
table  hs  Part  HI  of  Volume  1  of  these  Report*,  giving  v.  in  terms  of  M  and  4,.  One  feature  of  these  tables  and 
diagrams  should  be  noted  particularly:  namely,  the  rapid  rise  of  both  4/  *  and  Km  in  the  domain  of  the  "second” 
sShadoB*.  The  reader  may  recall  that,  for  every  emu-apex  angle  of  a  solid  cone,  two  values  of  A,/t  can  be  found 
which  correspond  to  the  same  velocity  of  undisturbed  stream  in  front  of  the  shock  wave.  Solutions  correspond* 
tag  to  the  huger  and  smaller  value  of  Ajc  leading  to  identical  velocities  of  free  stream  are  customarily  referred 
to,  respectively,  as  the  "first”  and  “second”  solutions  of  our  hydrodynamics!  problem.  As  long  as  we  limit 
ourselves  to  a  purely  axial  Sow,  both  solution*  are  equally  admissible  from  the  mathematical  point  of  view. 

1*.  Aa  aseaptioD  to  tide  are*  is  ostein  sohstio-ia  is  the  aanediai*  Mfefcborhood  ot  tlw  solid  case,  when  say  interpolation  tar  s 
heesms  tut  illfi  nf  Into* d  of  hwahteg  itows  oo  tofadrtiou  to  wysnsB^eps  —  which  would  han  boos  aatfasly  — portoon  for  *H 

nlw  II  dm  m  ItoMad  to  hridgs  Ik - tetooifcoad  of  tl - at  largo  iat - Such  cans  can  be  eaaily  ntcogaiaad  by  as  odd 

am  «f  than  In  as  the  hagtoato*  «f  tka  sospaetim  tobslstion. 

14.  Igft  to  F sera- Bay,  flu  tot.  Sne.,  S,  IT,  IM. 
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f  VYL*oi,  has revex,  arc  cotnc  lo  consider  the  first-order  effects  of  yew,  we  are  confronted  with  a  deferent  mostinp 

For  a  given  pair  of  t.  and  U /alt  the  value*  of  both  t/ 1  and  Kn  corresponding  to  a  "second"  motion  are  found 
to  be  aeveral  time*  as  Urge  aa  those  appropriate  for  the  "first”  one,  and  their  disparity  increases  rapidly  with 
increasing  Mach  number.  Yet  this  is  in  discord  with  our  initial  assumption  that  the  yawing  dr-ck  wave  Aoald 
remain  "dow”  to  its  ooo-yaw  position,  and  makes  it  therefore  questionable  to  what  extent  our  theory  is  actually 
applicable  to  the  "second"  solutions.  At  any  rate,  there  seems  but  tittle  room  for  doubt  tha  t  the  “second" 
i  solution*  for  a  supersonic  flow,  around  yawing  cones,  characterised  by  fairly  Urge  Mach  c embers  must  be 
J  unstable  to  a  high  degree.  Whether  or  not  this  has  anything  to  do  with  the  fact  that  shock  wavefr  corresponding 

'  to  these  “second”  solutions  have  never  ms  yet  been  observed  on  spark-range  photographs  (where  even  a  very 

■mall  yaw,  variable  in  time,  might  smooth  out  the  pressure  gradient  across  the  shock  and  th  is  considerably 
|  bhir  ita  features) ,  while  their  discoveries  have  from  time  to  time  been  reported  from  various  wni  t  tunnel*  (where 

£  the  model  can  be  maintained  in  a  fixed  position)  remains,  however,  to  be  clarified  by  future  investigations. 


DISCUSSION 

Before  proceeding  to  discuss  the  application*  of  the  foregoing  theory’  to  practical  cases,  a  recapitulation 
of  the  physical  assumptions  involved  in  it  may  be  in  order.  First,  we  have  tacitly  assumed  that  the  Mach 
number  (7/ai  of  free  stream  is  high  enough,  and  the  semi-apex  angle  6,  of  the  cone  is  small  enoi  -H,  for  the  non¬ 
yaw  case  to  give  rise  to  a  conical  shock  wave  attached  to  the  vertex  of  the  cone.  This  amounts  to  a  requirement 
that  U/ai  be  somewhat  larger  than  unity,  and  6,  less  than  56°,  the  exact  limits  being  given  in  Volume  1  of  these 
Reports.  Secondly,  the  angle  of  yaw  *  is  assumed  to  be  small  enough  for  quantities  of  the  order  ?tf  «*  and  higher 
powers  to  be  negligible,  and  the  shock  wave  is  assumed  to  be  "close"  to  ita  non-yaw  position  it  is  expected 
to  devUte  from  the  Utter  by  amounts  of  the  order  of  «).  Third,  the  how  reUtive  to  the  cone  is  assumed  to  be 
steady,  and  its  properties  to  be  independent  of  r  and  constant  over  co-axial  cones  (not  necessarily  circular) 
passing  through  the  same  vertex.  Lastly,  the  flow  is  supposed  to  be  adiabatic  along  each  streamline,  though 
not  iaentropie  over  the  whole  field.  That  is,  the  entropy  should  be  constant  —  except  for  a  discontinuous 
increase  at  the  shock  wave  —  along  each  streamline,  but  will  in  general  vary  from  One  streamline  to  another. 
It  can  be  show  n  that  this  amounts  to  a  neglect  of  the  effects  of  viscosity  and  heat  conduction  in  the  Bow  between 
the  shock  wave  and  the  cone.  The  air  will  then  behave  like  a  perfect  fluid  and  an  ideal  gas;  and  the  standard 
Rankine-Hugoniol  conditions  should  be  obeyed  at  the  shock  wave.  It  should,  however,  be  emphasised  that  at 
no  stage  of  our  analysis  did  we  assume  the  flow  to  be  irrotational. 

The  preceding  theory  is  concerned  with  an  idealised  projectile  consisting  of  an  infinite  cone  which  moves 
with  a  constant  speed  wed  yaw.  Consider,  instead,  an  actual  projectile  having  a  conical  head.  Its  speed  will  be 
changing  slowly,  and  its  yaw  may  vary  considerably ,  though  the  change  in  yaw  is  ordinarily  slow  when  compared 
with  the  velocity  of  transklian  of  the  projectile.  Under  such  circumstances,  th*  flow  of  air  around  it  can 
indeed  be  regarded  as  steady  to  a  high  degree  cf  approximation,  and  to  be  very  nearly  the  same  as  the  one  defined 
by  our  theory  —  provided  only  that  it  remains  locally  supersonic  in  the  whole  region  between  the  cone  and  its 
shock  wave.1*  Thus  for  an  actual,  fairly  stable  conical-headed  projectile  moving  at  a  requisite  supersonic  speed 
and  subject  to  a  small  yaw  «,  the  shock  wave  should  remain  a  circular  cone  of  the  same  semi-apex  angle  fw 
as  if  there  were  no  yaw,  but  yawed  at  an  angle  i,  proportions!  to  t  as  given  by  equation  (48-1)  with  respect 
to  the  direction  of  free  stream 

This  provides  an  experimental  .rack  on  the  theory.  By  means  of  spark  photography,  shadowgrams  can 
be  obtained  which  give,  in  effect,  the  projections  of  a  missile  in  supersonic  flight  and  of  its  accompanying  shock 
wave  upon  the  planes  of  strategically  oriented  photographic  plates.  Appropriate  measurements  of  a  pair  of 
such  plates  can  then  disclose  the  instantaneous  plane  and  magnitude  of  the  yaw.  Careful  investigations  of 
such  photographs,  taken  at  the  spark  range  of  Free  Flight  Aerodynamics  Branch  of  the  Ballistic  Research 
Laboratories,  U.S.  Army  Proving  Ground  at  Aberdeen,  Maryland,  of  a  15°-cone  characterised  by  the  Mach 

IS.  Without  thi»  wwrvstton.  the  affect  o f  the  truncation  of  the  cone  might  be  propagated  downstream  and  thua  modify  the  Sow 
scorned  the  eoeueal  head. 
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wliw  1.90  by  I.  t.  Segal  diadoeed  the  agreement  between  theory  and  observation*  to  be  very  satu- 

feetory  lor  values  of  «  up  to  -v  least  5*.  A  mors  extensive  series  of  observations  aiming  at  an  observational 
verification  of  Stone's  theory  iu  recently  undertaken  by  Professor  R.  W.  Lsdenburn  and  bis  collaborators  at 
the  Palmer  Physical  Laboratory  of  Princeton  University.  The  Princeton  investigators  used  a  supersonic  jet  of 
compressed  air,  characterised  by  a  Mach  number  of  1.70,  and  experimented  with  three  different  cones  of  semi- 
apex  ingin  of  10*,  *0*,  and  30s,  exposed  to  the  supersonic  stream  at  vsrious  angles  erf  mrimatam  (see  the 
Ftantaspiece).  la  a  memorandum  report,  kindly  communicated  to  the  writer  in  advance  of  publication.  Laden 
burg  and  Bemhader  concluded  that  .  .  "Qualitatively,  the  experimental  results  confirm  the  theoretical  expec¬ 

tation,"  while  "The  quantitative  agrnrmwt  is  .  .  .  within  the  relatively  large  uncertainties  of  the  measure¬ 
ment*.”  On  the  strength  of  this  observational  evidence,  from  both  the  spark  range  and  the  laboratory,  ve 
eastern!  that  Stone's  bob  rtneer  theory  of  a  Rpmook  flow  around  slightly  yawing  cooes  appears  to  provide 
quite  a  dose  description  of  the  respective  physical  phenomena  —  certainly  the  best  one  available:  so  tar  —  and 
the  numerical  data  presented  in  this  volume  should  enable  one  to  compare  its  consequences  with  the  observa¬ 
tions  for  any  particular  case  with  a  minimum  of  difficulty. 

A  knowledge  of  these  data  puts  us,  incidentally,  in  a  position  to  determine  quantitatively  the  error  of 
earlier  repreeeatation*  of  the  superaooic  flow  around  yawing  cones  based  on  linearised  equations,  of  which  the 
osa  by  Tsssn*  appear*  to  be  the  simplest.  Tsien,  in  effect,  extended  the  von  Kkrmhn-Moore's  linearised  approx¬ 
imation  to  non -axial  flow  around  slender  cones  to  the  order  of  accuracy  to  which  squares  and  higher  powers  of  < 
can  be  ignored.  The  detail*  of  his  investigation  cannot  be  reproduced  here;  but  in  what  follows  we  shall  briefly 
wvp™  <*«  *onwqwv» »*«  with  ♦*»»  twtonm*  «*  non-linear  th— wy  We  found  that,  according  to  the  latter, 

the  first-order  yaw  terms  should  cause  the  axis  of  the  shock  wave  to  depart  from  the  axis  of  the  solid,  and  give 
rise  to  a  force  acting  on  tire  con e  in  a  direction  perpendicular  to  its  motion  The  former  phenomenon  has  no 
counterpart  in  Tsien’*  theory,  since  the  linearised  equations  do  not  account  for  the  existence  of  any  shock  wave. 
For  the  normal-force  coefficient  acting  on  a  slightly  yawing  slender  cone  Tsien  derived,  however,  an  expression 
of  the  form 


/_  ootfi 


The  following  tabulation  shows  a  comparison  of  Trim’s  values  of  Kg  as  computed  from  the  foregoing 
equation  (37)  with  those  evaluated  from  (84)  and  based  on  the  outcome  of  our  numerical  integration*  of  Stone’s 
non  fiwcsr  equations  for  three  different  ernes  and  a  fairly  wide  range  of  Mach  numbers: 


a* 

t7<*« 

Kk  (Tsien) 

6,-10* 

Km  (Stone) 

Error 

0,4 

1.000 

0.771 

0.715 

-7.8% 

0,5 

1.408 

0.781 

0.710 

-1.7% 

0.0 

1B17 

0.084 

0.715 

4.8% 

0.7 

t-387 

0.048 

0.718 

1*B% 

0B 

SBOt 

0.545 

0.718 

•4.1% 

0B 

5.44a: 

0.400 

o.tos 

444% 

0.06 

9.100 

(imaginary) 

0.7M 

— 
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04 

1.188 

0731 

0.668 

-9.4% 

0.46 

1.338 

0.689 

0.673 

-3.8% 

0.5 

1.514 

0.665 

0.675 

1.5% 

0.6 

1.954 

0.593 

0,678 

18.7% 

07 

8.579 

0.511 

0.687 

85.6% 

C.S 

3.834 

0,410 

0.701 

4!  5% 

0.95 

16.855 

(imaginary) 

0  744 

.... 

5,-80® 


0.35 

1.817 

0.06* 

0.615 

-18.8% 

0.44 

1.467 

0.681 

0.689 

1.3% 

0.55 

1.871 

0.536 

0.040 

!«.*% 

0.65 

8.443 

0.448 

0.058 

31.3% 

0.75 

3.S69 

(imaginary) 

0.607 

- 

0.85 

5.546 

(imaginary) 

0.686 

- 

A  glance  at  this  tabulation  reveals  that  the  linearised  theory  would,  on  the  whole,  be  a  very  poor  substitute  for 
the  more  correct  non-linear  one.  In  particular,  according  to  the  linearised  theory  the  coefficient  A'*  should 
decrease  with  the  increasing  Mach  number  and  vanish  when  Ufa  i  =  cos  8„  whereas  according  to  the  non-linear 
theory  A*  behaves  in  exactly  the  opposite  way  and  would  not  vanish  even  should  the  Mach  number  become 
infinite.  The  departures  from  linearity  —  or,  in  physical  language,  from  isentropy  and  irrotationality  —  thus 
tarn  out  to  exert  a  significant  influence  on  the  supersonic  Sow  around  yawing  cones,  and  this  influence  is  found 
to  increase  very  rapidly  with  increasing  velocity  of  the  free  stream. 
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These  tables  give,  successively,  the  contributions  x,  y,  z  to  the  velocity  components  u,  r,  u-  taken  in  the 
direction  of  increasing  spheric*!  polar  coordinates  r.  9,  $  due  to  a  small  yaw,  as  functions  of  the  angular  variable 
9;  and  the  corresponding  proportional  changes  in  pressure  (a  f)  and  density  (£  p)  due  to  the  same  cause. 

The  velocities  are  expressed  in  terms  of  the  velocity  which  the  air  in  front  of  the  shock  wave  would 
attain  if  allowed  to  expand  adiabaticaily  into  s.  vacuum. 

A  sign  of  (ri)  following  the  value  of  the  Mach  number  .!/  =  ('  <ii  in  the  headings  indicate  a  'second' 
solution. 
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SURVEY  OF  THE  RESULTS 


The  individual  eolumnn  indicate: 

it,  =*  radial  velocity -component  ct  the  axial  flow  along  the  solid  surface; 
semi-apex  angle  of  the  shock  wave; 

M  =  Mach  number  (t.e.,  stream  velocity  divided  by  the  velocity  of  sound  in  the  undisturbed  air  in  front 
of  the  shock  wave); 

d  ^constant  specifying  the  change  of  entropy  across  streamlines; 

S/t  =rutio  of  the  yaw  S  of  the  shock  wave  to  the  yaw  *  of  the  solid  cone; 

Kd~ coefficient  of  head  drag  on  the  cone; 

X"jv  =  coefficient  of  normal  drag  on  the  cone. 
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DIAGRAMS  ILLUSTRATING  SOME  PHYSICAL  PROPERTIES 
OF  SUPERSONIC  FLOW  OF  AIR  AROUND  SLIGHTLY  YAWING  CONKS,  SUMMARIZED 

IN  TABLES  II,  III,  AND  IV 


Dtp  individual  diagrams  indicate: 

l.  Shock-to-Cone  Yaw  Ratios,  S/t,  plotted  against 
4.  Coefficient  of  Normal  Drag,  A'y.  plotted  against  u,. 

Different  types  of  broken  lines  in  the  diagrams  indicate,  respectively: 

1.  The  parabolic  limit,  i.e.  the  values  of  ti,  corresponding  to  the  minimum  Mach  numbers  for  which 
the  stream  between  the  cone  and  the  shock  wave  remains  entirely  supersonic; 

2.  The  elliptic  limit,  i.e.  the  values  of  u,  corresponding  to  the  maximum  Mach  numbers  for  which  the 
stream  between  the  cone  and  the  shock  wave  is  entirely  subsonic; 

:t.  The  separation  between  “first”  and  “second”  waves,  i.e.  the  velars  of  v,  corresponding  to  the  mini¬ 
mum  Mach  lumbers  for  which  a  conical  regime  of  flow  is  possible. 

4.  Infinite  Mach  number,  i,e.  the  values  of  «,  corresponding  to  Mach  number  infinity. 


SHOCK -TO- CONE  YAW  RATIOS 
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